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ON TWO DIMENSIONAL REGULAR LATTICES 
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Abstract. Darboux transformations for linear operators on regular two dimensional lattices are re- 
viewed. The six point scheme is considered as the master linear problem, whose various specifications, 
reductions, and their sublattice combinations lead to other linear operators together with the corre- 
sponding Darboux transformations. The second part of the review deals with multidimensional aspects 
of (basic reductions of) the four point scheme, as well as the three point scheme. 



1. Introduction 

Darboux transformations are a well known tool in the theory of integrable systems [72] [91] [48] . The 
classical Darboux transformation [15j deals with a Sturm-Liouvillc problem (the one dimensional sta- 
tionary Schrodingcr equation) generating, at the same time, new potentials and new wave functions 
from given ones thus providing solutions to the Korteweg-de Vries hierarchy of 1 + 1 dimensional in- 
tegrable systems. However, as it is clearly stated in [H], it was an earlier work of Moutard [71] which 
inspired Darboux. In that work one can find the proper transformation that applies to 2 + 1 dimensional 
integrable systems. Note that the initial area of applications of the Darboux transformations, which 
preceded the theory of integrable systems, was the theory of conjugate and asymptotic nets where the 
large body of results on Darboux transformations was formulated [HH EH EE 11041 [Ml E] ■ 

Most of the techniques that allow us to find solutions of integrable non-linear differential equations 
have been successfully applied to difference equations. These include mutually interrelated methods 
such as bilinearization method [49] and the Sato approach [17j [19] , direct linearization method [81] [80] , 
inverse scattering method Q] , the nonlocal 8 dressing method [TO] , the algebro-geometric techniques [53] ■ 

Also the method of the Darboux transformation has been successfully applied to the discrete integrable 
systems. The present paper aims to review application of the Darboux transformation technique for the 
equations that can be regarded as discretizations of second order linear differential equations in two 
dimensions (|3.ip . their distinguished subclasses and systems of such equations. 

While presenting the results we are trying to keep the relation with continuous case. However, we 
are aware of weak points of this way of exposition. The theory of discrete integrable systems |1031 |4"S] 
is reacher but also, in a sense, simpler then the corresponding theory of integrable partial differential 
equations. In the course of a limiting procedure, which gives differential systems from the discrete ones, 
various symmetries and relations between different discrete systems are lost. The classical example is 
provided (see, for example pT]) by the hierarchy of the Kadomsev-Pctviashvilii (KP) equations, which 
can be obtained from a single Hirota-Miwa equation — the opposite way, from differential to discrete, 
involves all equations of the hierarchy [73] . 

The structure of the paper is as follows. In section [5] we expose main ideas in the theory of Darboux 
transformations in multidimension. In Scction[3]we present the construction of the Darboux transforma- 
tion for the discrete second order linear problem — the 6-point scheme (|3.2j) — which can be considered 
as a discretization of the general second order linear partial differential equation in two variables (13.11) . 
Then we discuss various specifications and reductions of the 6-point scheme. We separately present in 
Section [3.41 the Darboux transformations for discrete self-adjoint two dimensional linear systems on the 
square, triangular and the honeycomb lattices, and their relation to the discrete Moutard transformation. 
Section [4] is devoted to detailed presentation of the Darboux transformations for systems of the 4-point 
linear problems, their various specifications and the corresponding permutability theorems. Section [5] 
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is dedicated to reductions of the fundamental transformation compatible with additional restrictions on 
the form of the four point scheme. Finally, in Section El we review the Darboux transformations for the 
3-point linear problem, and the corresponding celebrated Hirota's discrete KP nonlinear system. 

2. Jonas fundamental transformation and its basic reductions 

The main idea included in the paper has its origin in Jonas paper |51| where fundamental transforma- 
tion for conjugate nets has been presented. Neglecting the geometrical context of Jonas paper we only 
would like to say that the fundamental transformation acts on solutions of compatible system of second 
order linear differential equations on function ip : W™ — ► IR m 

^,x*xi = ■A >3 if) >x i +A Jz ip, x i , i,j = l,...,n i^j 

where coefficients A 13 and A 31 are functions of and for n > 2 they have to obey a nonlinear 

differential equation (compatibility conditions); we use notation that subscript preceded by coma denotes 
partial differentiation with respect to indicated variables. 

Every fundamental transformation ip i— > ip can be presented as composition of suitably chosen 



(1) radial transformation ip t— ► ip r 

(2) Combescure transformation ip T 



(3) radial transformation ip r <— > ip 

i' = ?VV 

We refer interested reader to Eisenhart book [30] for further details. As the practice shows the fun- 
damental transformation is indeed fundamental - in the sense that most (if not all) of the Darboux 
transformations are reductions of the fundamental transformation. 

Here we follow this basic idea discovered almost one hundred years ago. We start from presenting two 
dimensional difference operators L (and corresponding difference equations = 0, where \I/ is function 
of discrete variables m and n) together with the transformation (>?, L) i— ► (^, L) which is composition of 

(1) gauge transformation (^,i) i— > (^ r ,L r ) (sec subsection 13. 2j) 

L r = $oIo 6, * r = — 

where is a solution of L0 = whereas $ is a solution of L^<£> = where U denotes operator 
formally adjoint to operator L. We emphasize that this particular choice of functions and 4> 
is essential, for it guarantees existence of function \Ev in the next transformation 

(2) transformations (^> ri L r ) ^ (*f> r ,L r ) that take either the form (details are given in the text 
below) 



(2.1) 



or the form 



(2.2) 
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"Ai* r " 
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A 2 * r . 
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A_ 2 * r 



where Ai, A 2 denote forward difference operators Ai/(m, n) := /(m+1, n)—f(m, n), A 2 /(m, n) : 
/(to, n + 1) — /(to, n), whereas A_i, A_ 2 denote backward difference operators A_i/(m, n) := 
/(to - 1, n) - /(to, n), A_ 2 /(m, n)j= /(to, n - 1) - f(m, n). 
(3) gauge transformation (if? r ,L r ) 1— ► (^ r ,L r ) 



L r = Ro L o S, 



*' = 8 



where in general R and S are arbitrary functions (which should be specified when reductions or 
specifications of the transformations are considered) . 



6-point 



sect. 13.11 



• • • 

+B* (22 ) +2C*(i2) +G*(1) +H*( 2 ) = i 7 * 

aV'ixs +bip, vv +2cip, xy +gip, x +hip, y = ftp 



-4(i)%) +-43'(-i) +B (2 )*(2) +B*(- 2 )+ 

C(l)1 , (l-2)+C(2)*(-12) = -P* 

\(aii, :c +cip, v ), n . + {bi!, y +cip, n ), y = fip J 



Honeycomb • • 
See section 13.4.31 for details 



Specification C = \ 



Specification 
A = = B, C = 



'5-point 
self adjoint 



• • • sect. [3~4~2I 



4-point *• sect. [3~3llO 



*(12)+G*(i)+-ff*( 2 )=-F* 
1p,xy +gip,x +hlp, y = ftp 



Ribaucour reduction 



8(2)1'(2)+e*(-2)=J C '* 

Sublatticc approach 
on square lattice 



Sublatticc approach 
on quasi-regular 
rhombic lattice 



Specification 
F = 



See section 

\I/ lies on a quadric 
tp lies on a quadric 



discre te Moutard 
sect. 13X51 

*(12)+* = A/(*(l)+*(2)) 

tp,xy = flp 



adjoint discrete 
Moutard, sect. 15X21 



*(12)+* = Af( 2 )*(2 

1p,xy = ftp 



Moutard reduction 



discre te Goursa t 
sect. 13311 

'(2) , 



(12) = 
T (12) T 1(2) 



*(l)+*(2)- 



T (1) T (2) q 
2 T (1) T (2)- T (12) T 



Goursat reduction 



3-point 


sect. 1331161 




Kfftt (a) =0 


gip, x +hip, y = 





> 

3d 
C 

H 

> 
cc 

% 
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C 
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Figure 1. Interrelations between schemes discussed in the text. The figure include: names of the schemes used in the text, position of 
the points the equation relates, numbers of sections where the DT for the scheme is discussed, the equation and its natural continuum 
limit (extensions of DT to multidimension are known only for the equations in rectangle frames) 
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We pay special attention to subclasses of operators L that admit Darboux transformations. A trivial, 
from the point of view of integrable systems, examples are subclasses obtained by fixing a gauge (see 
section f3.2[) . Very important classes of operators are obtained by imposing conditions on the mentioned 
above functions O and <f>. For example in the case of self-adjoint 7-point scheme we discuss in section 
13.4.11 the relation is $ = 9. We refer to such procedures as to reductions. The further examples of 
reductions are given in sections 13.3.11 13.3.21 and [5] We reserve separate name "specification" to the 
cases when the operator is reduced, whereas the transformations remain essentially unaltered (i.e. no 
constraint on functions 9 and $ is necessary). Examples of specifications are given in sections 13.31 I3.4."2l 
and 13.51 

With the notion of "specification" another important issue appears, since continuous counterparts of 
the specifications presented here can be viewed as a choice of particular gauge and independent variables. 
The interesting (because not well understood) aspects of discrete integrable systems come from the fact 
there is no theory how to change discrete independent variables so that not to destroy underlying 
integrable phenomena. Sublattice approach, widely used in theoretical physics, can be regarded, to 
some extent, as counterpart of change of independent variables. To some extent, because the only case 
studied in detail from the point of view of Darboux transformations is the case of self-adjoint equation 
(|3.29p . which in the discrete case consists of Moutard case (section f3.3.2j) . self-adjoint case (section l3.4j) 
and their mutual relations (c.f. [391 138] ). Interrelations between results presented here are summarized 
in Figure 1. 

Large part of the paper is dedicated to the specification in which the matrix in (|2.ip is diagonal. In 
this case one can consider systems of four point operators defined on lattices with arbitrary number of 
independent variables. The multidimensional lattices are extensively discussed in sections |31 and [SJ where 
compact elegant expressions for superpositions of fundamental transformations are presented either. 



3. TWO DIMENSIONAL SYSTEMS 

In this section we present discretizations of equation (|3.1I) and its subclasses covariant under Darboux 
transformations. 

3.1. General case. Out of the schemes that can serve as a discretization of the 2D equation 

(3.1) atp, mx +bip, yy +2cip, xy +gip, x +hip, y = ftp 
the following 6-point scheme 

(3.2) A* (u) + B* (22) + 2C* (12 ) + G* ( i) + £T* ( a) = F<$ 

deserve a special attention (coefficients A, B etc. and dependent variable *B are functions on Z 2 , subscript 
in brackets denote shift operators, f(m, ri)n) = /(m+1, n), f(m, «)( 2 ) = f{Tn, n+1), f{m, n)(ii) = f{m+ 
2, n), /(to, n)(22) = f( m i n + 2) and /(to, n)(i 2 ) = /(to + 1, n + 1)). The scheme admits decomposition 
[{a\Ti + a 2 T 2 + a^)(f3iT 1 + /? 2 T 2 + k) + y]ip = (where T\ and T 2 are forward shift operators in 
first and second direction respectively) and therefore its Laplace transformations can be constructed 
(see [6H [El I105[ EOl Ell HI ES] and Section 14.2.31 for the notion of the Laplace transformations of 
multidimensional linear operators). What more important from the point of view of this review the 
scheme is covariant under a fundamental Darboux transformation |76j (the transformation is often 
referred to as binary Darboux transformation in soliton literature). Indeed, 

Theorem 3.1. Given a non-vanishing solution 9 of \ 3.2\) 

(3.3) AQ (11) + B9 (22) + 2C6 (12) + G9 (1) + HQ (2) = FO 
and a non-vanishing solution $ of the equation adjoint to equation iS.Oty 

(3.4) A(_x_!) $(_!_!) + B(_2-2)$(-2-2) + 2C(_ 1 _ 2 )#(_l_a) + + ff(-2)$(-2) = F$ 

(negative integers in brackets in subscript denote backward shift e.g. /(m,n)(_i) = f(m — l,n) etc.) the 
existence of auxiliary function P is guaranteed 

A_i($e ( i 2) P) = £(_2)$(-2)0(2) + B *©(22) + C(_l)$(-1)6(2) + G$9 (12) + tf$9( 2 ) 

^■ 0) A_ 2 ($6 (12) P) = -(A(_ 1) $ ( _ 1) 9 (1 ) + A$9 (11) + C ( _ 2) $ ( _ 2) 6 (1) + G$9 (12) + G$9 (1) ). 
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A 2 



Then equation 13. 2\) can be rewritten as 
(3.6) 

which in turn guarantees the existence of functions suc/i £/iai 
(3.7) 



'Ai(*')" 




A 2 (*'). 





^(-2)*(-2)®(2) 



)J=° 






[M 


PI 




A 2 ( 


S)J 



(C(_ 2) +P(-2))$(-2)0(l) 

-A ( _i)$ ( _i)e (1) (P(_i) - )$(_!) e (2) 

Assuming that matrix in \3. 7| ) is invertible on the whole lattice, i.e. 

D := [(P(-i) - C , (_i ) )(P(_2)+C(_2)) + -4(-i)-B(_2)]S(i)0(2)"I , (-i) ( I>(-2) / everywhere and finally on intro- 
ducing ^ via 



(3.8) 



* = — 

S 



and taking the opportunity of multiplying resulting equation by non-vanishing function R, we arrive at 
the conclusion that the function ^ satisfies equation of the form \3. 2\) but with new coefficients 



A _ ggfll)gg(llj 



D M ' 



RS, 



(3.9) 



12) / 9(11)^(1-2) (C+-P)(l-2) , 

G = -RS (1) (®£HiiJ! 



9(22) ^(-12) (C — P)(_i 
-Dm 



2)(c+p) (1 _2)+9(h)$a , e(2)^(-i)(c-p)(-i)+e(i)$(-i)A ( _i) 



C(i) 



H = -RS {2) ( e<22)j>( - 12) 



(C-P)(_i2)+9(22)$B , 9 ( i)$ ( _2)(C + P)(-2)+9(2)$(-2)B(_ 2) 



The family of maps \& i— > given by ()3.7p - (|3.8|) we refer to as Darboux transformations of equation 

HO). 



3.2. Gauge equivalence. We say that two linear operators L and L' arc gauge equivalent if one can 
find functions say <& and 8 such that L'=$oLo9 (where o stands for composition of operators). The 
idea to consider equivalence classes of linear operators (|3. 1[) with respect to gauge rather than single 
operator itself, goes to Laplace and Darboux papers J5H [T5] . In the continuous case it reflects in the 
fact that one can confine himself e.g. to equations (the so called affine gauge) 

atp, xx +blp,yy +i>,xy +g4>,x +hlp, y = 

or to (in this case we would like to introduce the name basic gauge) 

aip, xx +bip, yy +2cip, xy +k iy ip, x -k, x ip, y = 

without loss of generality. 

The Darboux transformation can be viewed as transformation acting on equivalence classes (with 
respect to gauge) of equation (|3.2p (compare [76]). Therefore one can confine himself to particular 
elements of equivalence class. Commonly used choice is to confine oneself to the affine gauge i.e. to 
equations (|3.2p that obeys 

A + B + 2C + G + H- F = 0. 

If one puts S = const in (|3.8[) (this condition is not necessary) then the above constraint is preserved 
under the Darboux transformation. One can consider further specification of the gauge 

A + B + 2C + G + H-F = 0, A ( _ 12) + 5 (1 _ 2) +2C + G (2) + H (1) - F {12) = 

This choice of gauge we would like to refer as to basic gauge of equation (|3.2p . Note that if the equation 
(|3.2p is in a basic gauge its formal adjoint is in a basic gauge too. 
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3.3. Specification to 4-point scheme and its reductions. In the continuous case due to possibility 
of changing independent variables one can reduce, provided equation p.l|) is hyperbolic, to canonical 
form 

i>,xy +gif>,x +H',y = ftp- 

In the discrete case similar result can be obtained in a different way. Taking a glance at (|3.9p we can 
notice that coefficients A, B and F transform in a very simple manner. In particular, if any of these 
coefficients equals zero then its transform equals zero too. Let us stress that in this case if we do so, 
we do not impose any constraints on transformation data, transformations remains essentially the same. 
First we shall concentrate on the case when two out of three mentioned functions vanish. 
If we put 

(3.10) A = 0, B = 0, C 



2 



then from equation (|3.9|) 

and one can adjust the function R so that 



A = 0, B = 



so we arrive at the 4-point scheme 

(3.11) * (12) = a# (x) +/3tf (2 ) + 7*- 

and its fundamental transformation c.f. [T0| . [3T | I7T | l3"S"] . 

We observe that the form of equation p. lip is covariant under the gauge L i— ► j^Lg and to identify 
whether two equations are equivalent or not we use the invariants of the gauge 

(3.12) K =^i, n=^i 

7(2) 7(1) 

Two equations are equivalent if their corresponding invariants k and n are equal [75] . 

3.3.1. Goursat equation. In this subsection we discuss the discretization of class of equations 

Iptxy = —1p,y +P 2 i> 

p 

which is referred to as Goursat equation. The discrete counterpart of Goursat equation arose from the 
surveys on Egorov lattices [99] and symmetric lattices [32] and can be written in the form [75j 

XT/ 9 < 2 \t/ Xtf, r (12)T g(2) lT , 

(3.13) W(i2) = W(i) + W( 2 ) W 

q r (i) T (2) 9 

where functions q and r are related via 

(3.14) g2 = r ( i)T (2) - W 



T 2 
T (l) 



The gauge invariant characterization of the discrete Goursat equation is either 

2 (1 +«(1))(1 + K (2)) 

(1 + «)(1 + « (12) ) 

or 

2 ( 1 +«(l))( 1 +«(2)) 

(1) { ' (l + n)(l + n (12) ) 

The Goursat equation can be isolated from the others 4-point schemes in the similar way that Goursat 
did it over hundred years ago [451 175] i.e. as the equation such that one of its Laplace transformations 
maps solutions of the equation to solutions of the adjoint equation. Therefore in this case if obeys 

(Emu 

a q ^ 2) ft x ft T{12)T g(2) ft 
(3-loj B(i2) = fe»(i) + B( 2 ) c> 

q r (i) T (2) 9 
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then its Laplace transformation 
(3.16) 



1 T (2) 



Ai0 



(2) 



1(2) T (12) 

is a solution of equation adjoint to equation (|3.13|) [75] . Equation (|3. 16|) is the constraint we impose on 
transformation data (i.e. functions 9 and $) in the fundamental transformation (|3.7[) . (|3.9[) (we recall 
we have already put A and B equal to zero). In addition if obeys (|3.15[) then 



(3.17) 



Ai ( —e z 

- T (2) 



A 2 



- 2 — (Aie) 2 

9 r (i) 



and as a result there exists function $ such that 



(3.18) 



Aitf = — (AiO) 2 , A 2 t? = — e 2 



'(2) 



Now it can be shown that one can put 



1 , 



*e (12) (p--) = -0 (ia) 



(3.19) 

and the transformation (|3 . 7[) takes form 
[Ai( 



<j>© (12) (p + -) = 



r (2) VAn?A 2 ^ 




'(2) 



(2) 



(3.20) 



A 2 (, 



i5 



/ T(2) ygTgSjg 





T(2) x/^? 

which is the discrete version of Goursat transformation 

f = T'd. 






AS). 

The transformation rule for the field r is 



3.3.2. Moutard equation. In this subsection we discuss discretization Moutard equation and its (Moutard) 
transformation [71] 

1p,xy = fll>. 

Moutard equation is self-adjoint equation (which allowed us to impose reduction <!> = in the continuous 
analogue of transformation (|3.7|) c.f. 73\). The point is that opposite to the continuous case there is no 
appropriate self-adjoint 4-point scheme 

We do have the reduction of fundamental transformation that can be regarded as discrete counterpart 
of Moutard transformation. Namely, class of equations that can be written in the form 



(3.21) 



*(12) + * = + * (2) ) 



we refer nowadays to as discrete Moutard equation. It appeared in the context of intcgrable systems 
in [H] and then its Moutard transformations have been studied in detail in [52] ■ Gauge invariant 
characterization of the class of discrete Moutard equations is [75] 



7l( 2 )n = 

Let us trace this reduction on the level of fundamental transformation. Putting 2C 
H =: MF the equations (|3.2[) . (|3.3p and (|3.4p take respectively form 



-F and G 



(3.22) 



F(*(12) +*) = G(*(l) +*(2)) 



(3.23) 



P(0(12)+0)=G(0(l)+0(2)) 



(3.24) F ( _i_a)$(_i-a) + F<5> = G ( _i)$(-i) + G ( _ 2) $ ( _ 2) 

The crucial observation is: if the function satisfies equation (|3.23[) then the function <I> given by 

(3.25) $ = i(0 (1) +0 (2) ) 
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satisfies equation (|3.24|) [75] . If we put 

e (1) - e (2) 



2P 



"Ai(e*)" 




"ee (1) 





A 2 (e*)_ 







-ee (2) 





"A — " 




a4. 



then equations (|3 . 5|) will be automatically satisfied. If in addition we put in (13. 8|) S = 9 then Darboux 
transformation (JSUJ takes form (c.f. [55]) 

(3.26) 

and it serves as transformation $ i— > ^ that maps solutions of equation (|3.21|) into solutions of another 
discrete Moutard equation 

(3.27) *(i2)+* = M(* (1) +* (a) ) 1 M= (1) ™ M 

W(12) W 

We recall that 9 an arbitrary non-vanishing fixed solution of the equation (|3.21|) 

(3.28) e (12) + e = M(o (1) + e (a) ) 

Every equation that is gauge equivalent to (|3.24p we refer to as an adjoint discrete equation Moutard 
equation, its gauge invariant characterization is [75| 

«(X2)«(1) = ™(12)™(2)- 

Equation (|3.23|) can be regarded as potential version of equation (|3.24[) 

3.4. Self-adjoint case. In the continuous case there is direct reduction of the fundamental transforma- 
tion for equation (|3.1[) to the Moutard type transformation for self-adjoint equation [75] 

(3.29) (aip )X +cip, y ), x +{bi/), y +ctp, x ), y = ftp. 

In this Section we present difference analogues of equation (|3.29p which allow for the Darboux transforma- 
tion. 9pposite to the continuous case the transformation will not be direct reduction of the fundamental 
transformation for the 6-point scheme (|3.2p . The self-adjoint discrete operators studied below are how- 
ever intimately related to discrete Moutard equation which provides the link between their Darboux 
transformations and the fundamental transformation for equation 



3.4.1. 7-point self-adjoint scheme. The following 7-point linear system 

(3.30) Al)*(l) + + #(2)*(2) + S* ( _2) + C(l)#(l-2) + C(2)*(-12) = F*, 

allows for the Darboux transformation [78j . 

Theorem 3.2. Given scalar solution of the linear equation (|3.30[) . then "J given as solution of the 
following system 

(Ar(§0)\ ( C0 ( _ 1} f _ 2) -e^Be + CG^)} /A_r| 
[6 - 6 > \A 2 ^0)J {0 { -d {A0+C0 { _ 2) ) -C0 ( _ 1) ©(_2) J vA_ 2 | 

satisfies the 7-point scheme (|3.30|) with the new fields given by 

(3 32) A °°W A n. °°^ B n C(-i-2)0(-i)0(-2) 




(3.33) F -- 

where 

(3.34) V = 0AB + { - 1) CA + 0(-2)CB. 

As it was shown in [35] the self-adjoint 7-point scheme (|3.30p can be obtained from the system 
of Moutard equations imposed consistently on quadrilaterals of the bipartite quasi-regular rhombic 
tiling (sec Fig. [2]), which is a particular case of the approach considered in [9]. Then the Moutard 
transformations can be also restricted to the triangular sublattice leading to Theorem [ 
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Figure 2. The bipartite quasi-regular rhombic tiling and the triangular and honeycomb lattices 

3.4.2. Specification to 5-point scheme. The 7-point scheme admits specification C — (alternatively one 
can put A = or B = 0) and as result we obtain specification to 5-point self-adjoint scheme [75]. 



(3.35) 



.4(1)*(1) + + B(2)*(2) + £*(-2) = J 7 *. 



The self-adjoint 5-point scheme and its Darboux transformation can be also obtained form the Moutard 
equation on the (bipartite) square lattice [39j . 



3.4.3. The honeycomb lattice. It is well known that the triangular and honeycomb grids are dual to each 
other (see Fig. [5]). Restriction of the system of the Moutard equations on the rhombic tiling to the 
honeycomb sublattice gives |38| the following linear system 



(3.36) 
(3.37) 



*-) 



1 



-*--*++- 

A "(1-2) 



C(12) 
1 



(* + 



(2) 



0, 



- f + ) = 0. 



Remark. Because "P - and ^ + satisfy separately the self-adjoint 7-point schemes p.30p . but with different 
coefficients, then the linear problem (|3.36p - (|3.37[) can be considered as a relation between two equations 
(|3.30p . This is the Laplace transformation between self-adjoint 7-point schemes studied in [551 1ST] . 

The corresponding restriction of the Moutard transformation gives the Darboux transformation for 
the honeycomb linear problem. 



Theorem 3.3. Given scalar solution (0 + ,0 ) of the honeycomb linear system (|3.36ll - p.37p then the 
solution ( v f + , \& _ ) of the system 

(3.38) 

(3.39) 

(3.40) 

with 1Z given by 
(3.41) 



= 4^[er_ 1 _ a) *7_ a) -er_ a) * ( _ 



(-12) 



K 

g(-2) 

K 



(-2), 



satisfies a new honeycomb linear system with the coefficients 



(3.42) 



Z-^tHe- - 



B = 



B, 



(-2) 

n 



(-1) (-1-2)' 



c ( i2) - ^o { _ 1} e { _ 2) . 
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3.5. Specification to 3-point scheme. We end the review of two-dimensional case with specification 
to the 3-point scheme i.e. to a discretization of first order differential equation 

gifj, x +hip, y = ftp. 

If we put 

A = 0, B = 0, F = 

then according to (13.91) 

A = 0, 5 = 0, F = 



It means that the fundamental transformation (|3.7p is also the Darboux transformation for the 3-point 
scheme O E2] 

2Ctf ( i2) + + # * (2) =0. 

This elementary scheme is the simplest one from the class considered here, but it deserves a special 
attention, because it leads to one of the most studied integrable discrete equation [35]. We confine 
ourselves to recalling briefly in Section [3 main results in this field. 
To the end let us rewrite the 3-point scheme in the basic gauge 

(3.43) (U( 2 ) - U(l))*(l2) + ("(l) - U)*(l) - (U(2) - w)*( 2 ) = 

4. The four point systems 

In this Section we present the Darboux transformations for the four point scheme (the discrete Laplace 
equation) from the point of view of systems of such equations, and the corresponding pcrmutability 
theorems. To keep the paper of reasonable size and in order to present the results from a simple 
algebraic perspective we do not discuss important relations of the subject to incidence and difference 
geometry |2JB ISS IZ1 US 1211 HS1 12Z1 QSl 1321 IS31 EM EHD [33 EU (see also [33 E3] and 
earlier works [§7J EE]), application of analytic [13 [33 [33 [32 [Ml [Ml [Ml [33 [33 H3 W\ and algebro- 
geometric J53IS3I11E2123I3I33I2311I1122] techniques of the integrable systems theory to construct 
large classes of solutions of the linear systems in question and solutions of the corresponding nonlinear 
discrete equations. 

In order to simplify discussion of the Darboux transformations for systems of the 4-point schemes we 
fix (without loss of generality [31]) the gauge to the affine one 

(4- 1 ) *(«) = + A m*w + ( x - - A m)*> 1 * i> 

where Aij : Z w — ► R are some functions constrained by the compatibility of the system (|4.ip . It is also 
convenient [10] to replace the second order linear system (|4.ip by a first order system as follows. The 
compatibility of (|4.1[) allows for definition of the potentials (the Lame coefficients) Hi : 1 N — > K such 
that 

(4-2) A «=^H7> 
The new wave functions xp i given by the decomposition 

(4.3) A,* = , 
satisfy the first order linear system 

(4.4) ^i = Qij(3)i>j, i^3, 

where the functions Qij , called the rotation coefficients, are calculated from the equation 

(4.5) A.//, O.JI,,. i^j, 

which is called adjoint to (|4.4p . Both the equations (|4. 5[) and (|4.4[) are compatible provided the fields 
Qij satisfy the discrete Darboux equations [TU] 

(4.6) A fc Q.y = Qik(k)Qkj , i^jj^kj^i. 



DARBOUX TRANSFORMATIONS 



11 



Corollary 4.1. Due to (|4.4|) the function xp i satisfies itself the four point equation 

(4.7) = + ^f^-^u) - ( jr^ (i - QmQm) ^ ^ j, 



w/izZe i/ie function satifies the adjoint of (|4.7|) m t/ie sense of (|3.4[) . 

The discrete Darboux equations imply existence of the potentials given as solutions of the compat- 
ible system 

(4.8) ean = i-Q m Q. m , 

Pi 

and yet another potential r such that 

(4.9) Pi = 2a. 

r 

In terms of the r-function and the functions 

(4.10) Tij = tQ.j , 

the meaning of which will be given in section 14.2.31 equations (|4.8|) and (|4.6p can be rewritten [37J 132 
in the bilinear form 

(4.11) t {ij) t = r (l) r b) - T ji(i) r ij{j) , i^j 

(4.12) Ty(*0 T = T- (fe )Tij + T ik ( k )T~kj , i^ j^k^i. 

4.1. The vectorial fundamental (bilinear Darboux) transformation. We start with a simple 
algebraic fact, whose consequences will be discussed throughout the remaining part of this Section. 

Theorem 4.2 ( [71\). Given the solution <p i : 7L N — > U, of the linear system (|4.4[) . and given the solution 
<p* : Z w — > W* , o/ £/ie adjoint linear system (|4.5p . These allow to construct the linear operator valued 
potential ft [</>,<£*] : Z N -> L(W,U), defined by 

(4.13) A;n[0,0*] = <&<g>0* (i) , i=l,...,N. 
7/W = U a?i<i i/ie potential Q is invertible, £l[<p, <p*] £ GL(W), then 

(4.14) fa : = n[0,0*]-Vi, i = i,...,Jvr ) 

(4.15) 0* = </)*n[0,0T\ i = l, 
satisfy the linear systems (|4.4p - (|4.5p correspondingly, with the fields 

(4.16) Qtf = Q«-(0i|n[0,0*] _1 l^i), ' v. i^j,. 

In addition, 

(4.17) O[0,0*] =C-fl[0,0*]- 1 , 
where C is a constant operator. 

Remark. Notice that because of f|4.3[) we have ^ = O [i/j, i?]. 

Corollary 4.3 ( |69j ) ■ The potentials pi and the r-function transform according to 

(4.18) ft = ft (i + ^ i) n[0,0*]-Vi). 

(4.19) f = rdetn[0,0*]. 

Applying the above transformation one can produce new compatible (afhnc) four point linear problems 
from the old ones. 

To obtain conventional transformation formulas consider [35] the following splitting of the vector 
space W of Theorem [ 



(4.20) W = E ® V © F , W* = 
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if 

(4.21) fc = I Ojj , <j>* =(0,0*,ip*) 
then, the corresponding potential matrix is of the form 

/Ie niW 

(4.22) si[4>,<t>*] = o n[o,e* 

\0 

and its inverse is 

/i E -n^rMfl,^]- 1 -n[</>,t/>*] + ri[?/>,0*]rj[0,0*]- 1 r2[0,'i/>* 

(4.23) n[0,0*] _1 = o n^,©*]- 1 -f}[0,0*]- 1 ^[0,'/'1 

\0 I F 

Let us consider the case of iT-dimensional transformation data space. V = M A , and F = R, E = R M , 
then the transformed solution * = Sl[ijj, H] of the four point scheme (recall that St = fi[i/>, i?]) up to a 
constant vector reads 

(4.24) * = * - n[v>, e*]n[e, 0*p 1 fi[0, h\, 

where the corresponding transformed solutions of the linear problem (|4.4|) and Hi of the adjoint 
linear problem (|4. 5|) are given by equations 

(4.25) ^^^-n^.^np.et 1 ^ 

(4.26) a< = f 4 - ejoiMT^P.^L 

and 

(4.27) Q y = Qij - e*n[e, e*]- 1 ^. 

The scalar (K = 1) fundamental transformation in the above form was given in [57j . 

Remark. To connect the above formalism to the results of Section [3] notice that given scalar solution 9i 
of the lineai' system (14.41) then the potential = O[0,H] is the scalar solution of the system (14.11) of 
second order linear equations. Moreover, given the solution 8*, of the adjoint linear system (|4. 5|) then 
the functions 

(4.28) a t = ^ 
satisfy [35] the linear system 

(4.29) Aj cTj = (Ay - 1) (cr j(j) - o- iW )) , i 7^ J. 
Equations (|4.29[) imply that the functions 

(4.30) $y = cr,^-) - <7j(ij), i ^ j, 
satisfy the corresponding equations 

(4.31) = A H®ij(-i) + A 3i$ij(-j) + (! - A ij(i) - A ji(j))$ij, i ^ j, 

adjoint of equations (|4.1[) . In the case M = 2 we obtain therefore the data of the transformation being 
the solution of the linear problem and the solution <E> = $12 of its adjoint, thus we recover results 
of Theorem 13.11 with the four point specification (|3.10p in the affine gauge; see [25J for more detailed 
description. Notice that in order to describe fundamental transformation in the second order formalism 
for M > 2 one should consider in addition the algebraic relations 

$*j = - 3>iii ^ij(-i-j) + ®jk(-j-k) = $ik(-*-fc)j 

which are consequences of definition (|4.30l) . 

It is important to notice that the vectorial fundamental transformation can be obtained as a super- 
position of K scalar transformations, which follows from the following observation. 
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Proposition 4.4 ([35]). Assume the following splitting of the data of the vectorial fundamental trans- 
formation 

(4-32) 6i= (%)' °* l= ( ^ °" ) ' 

associated with the partition W K = Wi Ka © M. Kb , which implies the following splitting of the potentials 



n[e a ,H] \ _ ( n[e a ,e* a ] n[e a ,e 



(4.33) n[e,H] ^ "^'~ J j , n[e,p] "^'^j n\e\ei} 

(4.34) n[v>,0*] = ( n[W a ] n[4>,ot] ). 

T/ien i/ie vectorial fundamental transformation is equivalent to the following superposition of vectorial 
fundamental transformations: 

1] Transformation \I> — ^^ a ^ with the data 9" , 6* ai and the corresponding potentials n[6 a ,H], Jl[0 a ,0*[, 

(4.35) * {a} = * - n[ij},e* a }ct[e a ,e* a \- 1 fi[e a ,H}, 

(4.36) ^ {a} = ^ - n[^, CM0 a , C] -1 ^, 

(4.37) H} a} = Hi - e* a n[o a , eT 1 ^, h]. 

2] Application on the result the vectorial fundamental transformation with the transformed data 

(4.38) e b{a} = e\ - np^jnpv:]- 1 *?, 

(4.39) e* b {a} = o% - e* ia n[o a , e* a ]- l n[o a , e* b \, 

and potentials 

(4.40) n[o\ h}^ = n[e b , h] - n[e b , e* a ]n[e a ,e* a ]- 1 si[e a , h] = n[e b{a} ,H^}, 

(4.41) n[e h , e* b }^ = n[e b , e* b } - n[8 b , d* a ]n[e a ,e* a ]- 1 n[e a ,e* b ] = n[e b{a \ e* b {a} ], 

(4.42) n[tf, 01}^ = nty, ei] - n[tf, CM0 a , ^[-^[r, 0:] = rt[^ a \ e* b {a} ], 



i.e. 



(4.43) * = * {a ' b} = * {a} - n[ip,8* b }^[n[d b ,8* b }^}- 1 n[d b ,H}^. 

Remark. The above formulas, apart from existence of the r-function, remain valid (eventually one needs 
the proper ordering of some factors) if, instead of the real field K, we consider [28 arbitrary division ring. 
Notice, that because the structure of the transformation formulas ([4.241) is a consequence of equation 
(|4.1?p then the formulas may be expressed in terms of quasi-determinants [35] (recall, that roughly 
speaking, a quasi-determinant is the inverse of an element of the inverse of a matrix with entries in a 
division ring). 

4.2. Reductions of the fundamental transformation. Let us list basic reductions of the (scalar) 
fundamental transformation. We follow the nomenclature of [35] which has origins in geometric termi- 
nology of transformations of conjugate nets [5T1 [131 [4TJI 11041 [Bl [B"3"I [4*T] . We provide also the terminology 
of modern theory of integrable systems O , where the fundamental transformation is called the 
binary Darboux transformation. All the transformations presented in this section can be derived [35] 
from the fundamental transformation through limiting procedures. 

4.2.1. The Levy (elementary Darboux) transformation. Given a scalar solution 9i of the linear problem 
the Levy transform of = fl[il>, H] is then given by 



(4.44) AW^-M^- 



14 



ADAM DOLIWA AND MACIEJ NIESZPORSKI 



The transformed Lame coefficients and new wave functions are of the form 

i-n[e t H\ , 

*'-%*>■ 

Within the nonlocal d-dressing method the elementary Darboux transformation was introduced in [TO] . 

4.2.2. The adjoint Levy (adjoint elementary Darboux) transformation. Given a scalar solution 9* of the 
adjoint linear problem (|4.5p the adjoint Levy transform £*(\&) of \& is given by 

(4.49) £'(*) = ¥-§n[^,0*]. 

The new Lame coefficients and the wave functions are of the form 



(4.45) 




(4 46) 




(4.47) 


CM 


(4.48) 


City,) 



(4.50) 




i 


(4.51) 




9* 


(4.52) 






(4.53) 







As it was shown in [35] . the scalar fundamental transformation can be obtained as superposition of 
the Levy transformation and its adjoint. The closed formulae for iterations of the Levy transformations 
in terms of Casorati determinants, and analogous result for the adjoint Levy transformation, was given 
in [67]. 

Remark. The description of the Levy transformation and its adjoint in the homogeneous formalism in 
the case of M = 2 is given in [35] . 

Remark. From analytic [lOj and geometric |35j point of view one can distinguish also the so called 
Combescure transformation, whose algebraic description is however very simple (the wave functions ip i 
are invariant). The Combescure transformation supplemented by the projective (or radial transforma- 
tion, whose algebraic description is also trivial [35|). generate the fundamental transformation. See also 
section [2] and 13.11 for generalization of the Combescure transformation. 

4.2.3. The Laplace (Schlesinger) transformation. The following transformations does not involve any 
functional parameters, and can be considered as further degeneration of the Levy (or its adjoint) reduc- 
tion. The Laplace transformation of \& is given by 



(4.54) £. iW = *_^. ) i + 

The Lame coefficients of the transformed linear problems read 
(4-55) C ij (H i ) = ^, 



(4.56) C ij (H j )=(Q ij A 



El 

Qij 



(-n 



(4.57) dj (H k ) = H k - Q^Hj , k?i, j, 
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and the new wave functions read 

(4-58) £ i .( V ,.) = _A i V, i + M^T/, i , 

(4-59) dj (</>,•) = --^ , 



The Laplace transformations satisfy generically the following identities 



(4.60) £ij&k)=i>k-7r*l>i, 



(4.61) 


dj < 




= id, 


(4.62) 


Cjk < 




= -Cjfe 


(4.63) 


£fci < 


3 C i:i -- 


= £fcj 



The Laplace transformation for the four point affine scheme was introduced in j^D] following the geometric 
ideas of [98j and independently in [88| using the factorization approach. The generalization for systems 
of four point schemes (quadrilateral lattices) was given in |35j . 

Remark. As it was shown in [51] the functions r»j defined by equation (|4.10p are r-functions of the 
transformed four point schemes 

(4-64) 7ii=Aj(r), 
which, due to (]4. 6 1 [) . leads [5Dj to the discrete Toda system [15] . 

5. Distinguished reductions of the four point scheme 

In this section we study (systems of) four point linear equations subject to additional constraints, 
and we provide corresponding reductions of the fundamental transformation. The basic algebraic idea 
behind such reduced transformations lies in a relationships between solutions of the linear problem and 
its adjoint, which should be preserved by the fundamental transformation (see, for example, application 
of this technique in |86| to reductions of the binary Darboux transformation for the Toda system) . Some 
results presented here have been partially covered in Section [3.31 but in a different setting. 

5.1. The Moutard (discrete BKP) reduction. Consider the system of discrete Moutard equations 
(the discrete BKP linear problem [321 [85]) 

(5.1) * W) -* = ^(* W -* W ) 1 l<i<j<N, 

for suitable functions : Z w — > R. Compatibility of the system implies existence of the potential 
t b : Z N — * K, in terms of which the functions can be written as 

(5-2) F ij = ^l, 

which satisfies system of Miwa's discrete BKP equations [73] 

The discrete Moutard system can be given [25] the first order formulation (|4.4p - (|4.5|) upon introducing 
the Lame coefficients 

t b 

(5.4) H i = (_i)£ fe<i «HJz^ 
and the rotation coefficients (below we assume i < j) 

(5.5) Qa = -(_1)£«»<," 

(5.6) Qji = -(-l)^k<j n 



T (i~3) 




T (-j) 


T (-j) 




T B 


T {-m 




\ T (~i) 


T B 




1 T B 
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which in view of (|4.8|) . gives the familiar relation between the r- functions of the KP and BKP hierarchies 
(5.7) t=(t b )\ 

The corresponding reduction of the fundamental transformation was given in |26] , where also a link with 
earlier work [55] on the discrete Moutard transformation has been established. 

Proposition 5.1 ([26|). Given solution 6i : Z w — > M. of the linear problem (|4.4|) corresponding 
to the Moutard linear system (|5.ip and its first order form (|5.4[) - (|5.6|) . Denote by ® = Q[6,H] the 
corresponding potential, which is also new vectorial solution of the linear problem (|5.1[) . 
1 ) Then 

-B 



(5.8) 9* = (_l)E fe<i - fe _y) (0 T t } 



provides a vectorial solution of the adjoint linear problem, and the corresponding potential £l[6, 6*] allows 
for the following constraint 

(5.9) n[e,e*] + n[e,e*f = 2©® © T 

2) The fundamental vectorial transform \& o/\f , given by <|4.24[) with the potentials ft restricted as above 
satisfies Moutard linear system (|5.ip and can be considered as the superposition of K scalar reduced 
fundamental transforms. 

Remark. In the scalar case and for N = 2 inserting in equations (|4.28[) and (|4.29[) the functions Hi and 
9*i given by (|5.4p and (|5.8|) correspondingly, we obtain the formula (|3.25|) in the gauge F = 1 (modulo 
the corresponding change of sign). 

Notice that given then, because of the constraint (|5.9p . to construct £1(0,0*) we need only its 
antisymmetric part S^®!®), which satisfies the system 

(5.10) A j 5(e|e) = © w ®e t -e®e^ ) . 

This observation is the key element of the connection of the above reduction of the fundamental trans- 
formation with earlier results [55] on the vectorial Moutard transformation for the system (|5.ip . where 
the formulas using Pfafnans were obtained (recall that determinant of a skew-symmetric matrix is a 
square of Pfaffian). In particular, the transformation rule for the r B -function can be recovered 

(t b Pf 5(01©), K even 

(5 ' u) ffl= wf° ; e \y 

{ \ e s(0|e> ) 

5.2. The symmetric (discrete CKP) reduction. Consider the linear problem subject to the con- 
straint [35] which arose from studies on the Egorov lattices [SS] 

(5.12) PiQji(i) = PjQij(j) i i 7^ j ■ 

Then the discrete Darboux equations ()4.6|) can be rewritten |100j in the following quartic form 

, g ( TT (vk) -T{i)T{jk) -T { j)T(ik) -T"(A;) T (y)) 2 + 4 ( r (i) r (j) T (fc) r (yfe)+' rT fe) r (jfc)' r (ifc)) = 

A ( T (i) T (jk) T U) T m + m T Ukmk)T(ij) + T (j) T (*fc) T W r (y))' 

which can be identified with equation derived in [53j in connection with the star-triangle relation in 
the Ising model. According to [100J, the above equation can be obtained from the CKP hierarchy via 
successive application of the corresponding reduction of the binary Darboux transformations. 

Construction [69j of the reduction of the fundamental transformations which preserves the constraint 
(|5.12p makes use the following observation. 

Lemma 5.2 (|32|). The following conditions are equivalent: 

1) The functions Qij, Pi satisfy constraint (|5.12p .- 

2) Given a nontrivial solution <p* of the adjoint linear problem (|4.5p then 
(5-14) cb t = Pl (cj>* {l) ) T 
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provides a solution of the linear problem 
3) The corresponding potential Sl[<f>, </>*] allows for the constraint 

(5.15) ^[</>,01 T = n[4>,4>*]. 

Theorem 5.3 ( [69| ) ■ When the data 9i, 9*, £l[9,9*] of the fundamental transformation satisfy condi- 
tions (|5.14[) - (|5.15[) then the new functions Qij, pi found by equations (|4.27[) . (]4.18p are constrained by 
(l5~T2l) . i.e. 

PiQji(i) = PjQij(j) ) * 7^ i ■ 

The corresponding permutability principle has been proved in [27). 

Remark. To connect the symmetric reduction with the Goursat equation and the corresponding trans- 
formation (see Section [333} notice that because of Corollary 14.11 and equations (|4.8[) - (|4.9p the function 
if) 2 satisfies equation (|3.13[) with q = Q2i(i)- Then, due to the same Corollary and condition (|5.14D . 
given scalar solution 62 of the four point equation of ip 2 then the function 

( rAi0 2 \ 



,M?21(l)/ (2) \ T ( 2 ^21(1)/ (2) 

satisfies its adjoint; in the lat equation we have used the linear problem (|4.4p and equation (|4.9 



5.3. Quadratic reduction. Consider the system of four point equations (|4.1| which solution \& is 
subject to the following quadratic constraint 

(5.16) * T Q* + a T * + c = ; 

here Q is a non-degencrate symmetric matrix, a is a constant vector, c is a scalar. 

Remark. Notice that unlike in two previous reductions we fix (by giving the quadratic equation) the 
dimension of \&. 

Double discrete differentiation of equation (|5.16[) in i ^ j directions gives, after some algebra, the 
condition 

(5.17) xpf^QxPj+xpj^QxP^O, 

analogous to that obtained in [31] in order to characterize circular lattices [7] U3J . It implies [33] , in 
particular, that iftfQxpi satisfy the same equation (|4.8[) as the potentials pi. 

As in two above reductions, the quadratic condition allows for a relation between solutions of the 
linear system (|4.4p and its adjoint (|4.5p . The following proposition can be easily derived from analogous 
results of [21], where as the basic ingredient of the transformation was used the potential 9*}, but 
we present here its direct proof in the spirit of corresponding results found for the circular lattice [571168] . 

Proposition 5.4. Given a nontrivial solution 9* of the adjoint linear problem (|4.5p corresponding to 
the system of four point equations (|4.ip which solution is subject to the constraint (|5.16p then 

(5.18) Oi = m, 6*] (i) + 8*]) T Qip % 
provides a solution of the corresponding linear problem (|4.4p . 

Proof. After some algebra using the equations satisfied by Y* and xp i one gets 

(5.19) A O, - Qij{j)6j = fffcfi QH-, + ^J W QV, 
which vanishes due to (|5.17p . □ 

The following result gives the discrete Ribaucour reduction of the fundamental transformation. 

Proposition 5.5 ([21]). Given solution 9* : Z N — > M. K of the adjoint linear problem (|4.5p corresponding 
to the quadratic constraint (15.161) . 
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1) Then the potentials £l[ij),9*\, fl\0,H\ and tt[9,9*], where 9i is the solution of the linear problem 
(|4.4p constructed from 9* by means of formula (]5.18p . allow for the constraints 

(5.20) n[9,H} T = 2ip T Qn[tp 1 9*} + a T n[ip 1 9*} , 

(5.21) fl[9, 9*} + fl[9, 9*} T = 2ft[ip, 9*} T Q^[tp, 9*] . 

2) The Ribaucour reduction of the fundamental vectorial transform St of \I r , given by (|4.24[) with the 
potentials ft restricted as above satisfies the quadratic constraint (|5.16p and can be considered as the 
superposition of K scalar Ribaucour transformations. 

As it was explained in [21], the Ribaucour transformations [57] of the circular lattice [7j [13] can 
be derived from the above approach after the stereographic projection from the Mobius sphere. The 
superposition principle for the Ribaucour transformation of circular lattices was derived also in [68] . 

6. The three point scheme 

In this final Section we present the vectorial Darboux transformations for the three point scheme (|6.1[) . 
The corresponding nonlinear difference system (|6.4p . known as the Hirota-Miwa equation, is perhaps 
the most important and widely studied integrable discrete system. It was discovered by Hirota [49] , who 
called it the discrete analogue of the two dimensional Toda lattice (see also [EE]), as a culmination of his 
studies on the bilinear form of nonlinear integrable equations. General feature of Hirota's equation was 
uncovered by Miwa [73] who found a remarkable transformation which connects the equation to the KP 
hierarchy [T§]. The Hirota-Miwa equation, called also the discrete KP equation, can be encountered in 
various branches of theoretical physics (92] EU H09j and mathematics [102] 160] [56] . 

Consider the linear system [18] 

(6.1) * w - = Uy*, i^j<N, 

whose compatibility leads to the following parametrization of the field Uij in terms of the potentials r, 

(6.2) r i(j ) = nllij, 
and then to existence of the r function 

(6.3) r . = (_!)£»<* 

T 

and, finally to the the discrete KP system [13 EI] 

( 6 - 4 ) 4) T m - T ") T m + T (Vfe)< i < j < k. 

The same nonlinear systems arises from compatibility of 

(6-5) *^-*J i) =E^*^. )j i^j, 

called the adjoint of (|6.1[) . 

We present the Darboux transformation for the three point scheme in the way similar to that of 
Section |4~T1 following the approach of [S2], see however early works on the subject [M] 194]. 

Theorem 6.1. Given the solution $ : 7L N — * U, of the linear system (|6.ip . and given the solution 
<&* : 1 N — > W*, of the adjoint linear system (|6.5p . These allow to construct the linear operator valued 
potential f2[*, **] : Z N -> L(W,U), defined by 

(6.6) A^I*,**] =*®** i)) i = l,.,.,N. 
IfW = U and the potential fl is invertible, **] G GL(W), then 

(6.7) * = 

(6.8) ** = $*n[$,**]- 1 , 

satisfy the linear systems (]6.ip and (|6.5p , correspondingly, with the fields 

(6.9) da = Uj - (^♦nt*,**]- 1 *)^ + (**«[*, s*]- 1 *)^. 

In addition, 

(6.10) «[$,$*] =C-fl[*,$*]- 1 , 
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where C is a constant operator. 

The transformation rule for the potentials rj reads 

(6.11) fi = n(i - 

while using the technique of the bordered determinants [50j one can show that |82j 

(6.12) f H = t h detn[*,#*]. 

More standard transformation formulas arise, when one splits, like in Section [4~T1 the vector space W 
of Theorem 16.11 as follows 

(6.13) W = E®¥ffiF, W*=E*®V*0F* ; 
if 

/*\ 

(6.14) *=l©j, ** = (0, ©*,**) . 

Then the corresponding potential matrix and its inverse have the structure like those in Section 14. 1[ 
which gives [52] 

(6.15) * = * fl[*, ©">[©, 0*]- 1 ©, 

(6.16) = 0*f2[0,0*] ^[Q,**], 

(6.17) Uu = Un - (©*«[©, e*]- 1 ©)^) + (©*«[©, e*]- 1 ©)^. 

Notice that one can consider [TH [501 [H] the three point linear problem in associative algebras. Then 
the structure of the transformation formulas (|6.10|) implies the quasi-determinant interpretation |44j of 
the above equations. 

Remark. The permutability property for the Darboux transformations of the three point scheme |84| 
can be formulated exactly (cancel subscripts i) like Theorem 14.41 

Finally we remark that the binary Darboux transformation for the three point linear problem can 
be decomposed |82j into superposition of the elementary Darboux transformation and its adjoint, which 
can be described as follows. Given a scalar solution O of the linear problem (|6.1|) then the elementary 
Darboux transformation 

(6.18) Z>(*) = * w - 

(6.19) ©(*•) = in[e, *•] , 

(6.20) V(t) = 6r 

leaves equations (|6.1[) and (|6.5p invariant. 

Analogously, given a scalar solution O* of the linear problem (|6.5|) then the adjoint elementary 
Darboux 

(6.21) = — n[e*,*], 

(6.22) P*(**) = -**_. t) + ^-!l** ; 

(6.23) X>*(t) = 6*t 

leaves equations (|6.1j) and ()6.5p invariant. 

The above transformations allow for vectorial forms, which can be conveniently written [82( 190] in 
terms of Casorati determinants. 
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7. Summary and open problems 

In the paper we aimed to present results on Darboux transformations of linear operators on two di- 
mensional regular lattices. To put some order into the review we started from the six point scheme (and 
its Darboux transformation) as the master linear problem. The path between its various specifications 
and reductions has been visualized in Figure [TJ We considered also in more detail the corresponding 
theory of the Darboux transformations of the systems of the four point schemes and their reductions. 
Finally, we briefly discussed the multidimensional aspects of the three point schemes. It is worth to 
mention that for systems of the three or four point schemes the Darboux transformations can be inter- 
preted as a way to generate new dimensions. This is very much connected to the pcrmutability of the 
transformations, which is a core of integrability of the corresponding nonlinear systems. 

Separate issue touched here is such extention of four point schemes that can be regarded as an analogue 
of discretization of a diferential equation in arbitrary paramctrization. More precisely we discussed here 
such an extension of general (AKP) case (section 13. lj) and Moutard-sclfadjoint (BKP) case (sections 
13.3.21 and 13 ,4[) . Goursat and Ribaucour reductions have not been investigated from this point of view. 
The multidimensional schemes that mimic equations governing conjugate nets (and their reductions) 
with arbitrary change of independend variables have been not exploited either. 

We also mentioned another approach to the problem of construction of the unified theory of the 
Darboux transformations. It consists in isolating basic "bricks" in order to use their combinations to 
construct more involved linear operators together with their Darboux transformations. Such an idea has 
been applied to derive, starting from the Moutard reduction of the four point scheme, the self-adjoint 
operators on the square (the five point scheme), triangular (the seven point scheme) and the honeycomb 
grids. Recently it was shown in |30j that the theory of systems of four point linear equations and their 
Laplace transformations follows from the theory of the three point systems. This means, that also the 
transformations of the four point scheme can be, in principle, derived from the three point scheme. An 
open question is if the six point scheme and its transformations can be decomposed in a similar way. 

Finally, we would like to mention possibility of considering (hierarchies of) continuous deformations 
of the above lattice linear problems, which would lead to (hierarchies of) discrete-differential integrable 
- by construction — equations. Some aspects of deformations of the self-adjoint seven and five linear 
systems have been elaborated in [95l HH [96] . 
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